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TAGGED PARTICLE IN SYMMETRIC SIMPLE EXCLUSION 
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Abstract. We prove a nonequilibirum central limit theorem for the position 
of a tagged particle in the one-dimensional nearest-neighbor symmetric simple 
exclusion process under diffusive scaling starting from a Bernoulli product 
measure associated to a smooth profile po : R — » [0, 1]. 



1. Introduction 

The asymptotic behavior of a tagged particle appears as one of the central prob- 
lems in the theory of interacting particle systems and remains mostly unsolved. 

The first important result on the position of a tagged particle in the diffusive 
scaling is due to Kipnis and Varadhan [3]. By proving an invariance principle 
for additive functionals of reversible Markov processes, Kipnis and Varadhan de- 
duced an equilibrium central limit theorem for the position of a tagged particle in 
symmetric simple exclusion processes. This result was extended by Varadhan [TO] 
for mean-zero asymmetric exclusion processes, through an invariance principle for 
Markov processes with generator satisfying a sector condition; and by Sethuraman, 
Varadhan and Yau [5] to asymmetric exclusion processes in dimension d > 3, re- 
laxing the sector condition by a graded sector condition. In these three contexts 
the authors prove that 

X t N 2 — E[X tN 2] 
N 

converges in law, as N f oo, to a Brownian motion with diffusion coefficient given 
by a variational formula. Here X t stands for the position of the tagged particle at 
time t. 

The noncquilibrium picture is much less clear. Even a law of large numbers for 
a tagged particle starting from a Bernoulli product measure with slowly varying 
parameter seems still out of reach. Rezakhanlou |S] proved a propagation of chaos 
result which states that the average behavior of tagged particles is described by 
diffusion process. A large deviations from this diffusive limit in dimension d > 3 
was obtained by Quastel, Rezakhanlou and Varadhan 0] . 

We prove in this article the first nonequilibrium central limit theorem for a tagged 
particle. Consider the one-dimensional nearest neighbor symmetric situation. In 
this context, as already observed by Arratia pQ, the scaling changes dramatically 
since to displace the tagged particle from the origin to a site A > 0, all particles 
between the origin and A need to move to the right of A. This observation relates 
the asymptotic behavior of the tagged particle to the hydrodynamic behavior of 
the system. The correct scaling for the law of large numbers should therefore be 
X t N2 /N and we expect {X t ^2 — E[X t j^2\)JyN to converge to a Gaussian variable. 
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The central limit theorem in equilibrium was obtained by Rost and Vares [J for a 
slightly different model. They proved that for each fixed t > 0, X t ^i />/N converges 
to a fractional Brownian motion Wt with variance given by E[W 2 ] = at 1 / 2 . We 
extend their result to the noncquilibrium case. 

The idea of the proof is to relate the position of the tagged particle to the well 
known hydrodynamic behavior of the symmetric exclusion process. Since particles 
cannot jump over other particles, the position of the tagged particle is determined 
by the current over one bond and the density profile of particles. Therefore, a 
nonequilibrium central limit theorem for the position of the tagged particle follows 
from a joint central limit theorem for the current and the density profile. Since the 
current over a bond can itself, at least formally, be written as the difference between 
the mass at the right of the bond at time t and the mass at time 0, a central limit 
theorem for the position of the tagged particle should follow from a nonequilibrium 
central limit theorem for the density field. This is the content of the article. 

There are three main ingredients in the proof. In Section|3we present a nonequi- 
librium central limit theorem for the current over a bond and show how it relates to 
the fluctuations of the density field. In sectionEJwe obtain a formula which relates 
the position of the tagged particle to the current over one bond and the density 
field. Finally, in Section we present a sharp estimate on the difference of the 
solution of the hydrodynamic equation and the solution of a discretized version of 
the hydrodynamic equation. 

2. Notation and Results 

The nearest neighbor one-dimensional symmetric exclusion process is a Markov 
process on {0, 1} Z which can be described as follows. Particles are initially dis- 
tributed over Z in such a way that each site is occupied by at most one particle. A 
particle at a site x waits for an exponential time and then jumps to x ± 1 provided 
the site is vacant. Otherwise the jump is suppressed and the process starts again. 

The state space of this Markov process is denoted by X — {0, 1} Z and the 
configurations by the Greek letter 77, so that r)(x) = 1 if site x is occupied for the 
configuration 77 and otherwise. The generator Ljq of the process speeded up by 
iV 2 is given by 

(L N f)( v ) = N 2 Y^[f(^' X+1 v)-f(v)}, 

where a x,x+1 r) is the configuration obtained from 77 by interchanging the occupation 
variables r)(x) and rj{x + 1): 

!r](x + 1) if Z = x, 
r](x) ifz = x + l, 
77(2) otherwise. 

For each configuration 77, denote by 71(77) the positive measure on R obtained by 
assigning mass iV -1 to each particle: 

7r(7?) = N~ 1 ^2r](x)5 x / N 

and let 7r t = 7r(?7t). 



NONEQUILIBRIUM CENTRAL LIMIT THEOREM 



3 



Fix a profile po : R — * [0, 1] with the first four derivatives limited. Denote by 
"^(.) the product measure on X associated to po: 

< ( . 3 {r?,^) = l} = Po(x/N) 

for x in Z. For each N > 1 and each measure p on X , denote by P M the probability 
on the path space D(M. + ,X) induced by the measure p and the Markov process 
with generator Lm- Expectation with respect to P M is denoted by E M . Note that 
we omitted the dependence of the probability P p on N to keep notation simple. 
This convention is adopted below for several other quantities which also depend on 
N . The hydrodynamic behavior of the symmetric simple exclusion process is trivial 
and described by the heat equation. 

Theorem 2.1. Fix a profile po : R — > [0,1]. Then, for all time t > 0, under 

P„« the sequence of random measures Tr t converges in probability to the absolutely 

po() 

continuous measure p(t, u)du whose density p is the solution of the heat equation 
with initial condition po-' 

dtp = Ap 

P(0, ■)=/*(•)• ( ' 

Here and below, A stands for the Laplacian. 

This theorem establishes a law of large numbers for the empirical measure. To 
state the central limit theorem some notation is required. For k > 0, denote by TLk 
the Hilbert space induced by smooth rapidly decreasing functions and the scalar 
product < ■, ■ >k defined by 

<f,g>k = <f, (x 2 -A) k g>, 

where < ■, ■ > stands for the usual scalar product in M. d . Notice that TCq = L 2 (R d ) 
and denote by H-k the dual of 7ik- 

Let pf(x) = E„n [% . A trivial computation shows that pf(x) is the solu- 

po(-) 

tion of the discrete heat equation: 

dtP?{x) = A N p?(x) , . 
f$(x)=po(x/N), (Z - Z) 

where (A N h)(x) - N*Y, V ,\ V - X \=i[h(y) ~ h(x)}. 

Fix k > 4 and denote by {Y t N ,t > 0} the so called density field, a H_fc-valued 
process given by 

Y t N (G) = ^^G^V)Mz)-pf(z)} 

for G in Hk- Denote by Qn the probability measure on the path space D(R+, H-k) 
induced by the process Y t N and the measure ^(.)- Next result is due to Galves, 
Kipnis and Spohn in dimension 1 and to Ravishankar 5 in dimension d > 2. 

Theorem 2.2. The sequence Qn converges to Q, the probability measure con- 
centrated on C(R + ,7i_fc) corresponding to the Orsntein-Uhlenbeck process Y t with 
mean zero and covariance given by 

E[Y t (H)Y s (G)} = [(T t _ s H)G X s- f dr f (T t - r H) (T s _ r G) {d r \r — A\ r } 

JK JO JM 

for < s < t and G, H € TLk ■ In this formula, {T t : t > 0} stands for the semigroup 
associated to the Laplacian and \s for the function x{s,u) = p(s, u) [1 — p(s, u)] . 
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Note that in the case of the heat equation, d r \r — A%r = 2(d x p) 2 . Also, in the 
equilibrium case, \ lSi constant in space and time so that the second term vanishes 
and we recover the equilibrium covariances. Finally, integrating by parts twice the 
expression with A\ r , we rewrite the limiting covariances as 

E[Y t (H)Y s (G)) = ( (T t H)(T s G) X o + 2 f dr f {VT t _ r H) (VT s _ r G) Xr , (2.3) 

JR JO JR 

where V/ is the space derivative of /. 

We examine in this article nonequilibrium central limit theorems for the current 
through a bond and the position of a tagged particle. For a bond (x,x+l), denote 
by J x ,x+i(t) the current over this bond. This is the total number of jumps from 
site x to site x + 1 in the time interval [0, t] minus the total number of jumps from 
site x + 1 to site x in the same time interval. 

Theorem 2.3. Fix u in R and let 

where xn — [uN\. Then, for every k > 1 and every < t\ < ■ ■ ■ < tk, (Z^, . . . , Z^ k ) 
converges in law to a Gaussian vector (Z tl , . . . , Z tk ) with covariance given by 

E[Z s Z t ] = ( dx P[B S < x] P[B t < x] X o(x) 

J —oo 
P oo 

+ / dxP[B s >x]P[B t >x] X o(x) 

JO 

/•s />oo 

+ 2 dr dxp t - r (0,x)p s - r (0,x)xr(x) 

JO J-oo 

provided s < t and u = 0. In this formula, B t is a standard Brownian motion 
starting from the origin andp t (x,y) is the Gaussian kernel. 

By translation invariance, in the case u ^ 0, we just need to translate x by — u 
in the covariance. 

Let H = 1{[0, oo)}. The covariance appearing in the previous theorem is easy to 
understand. Formally the current N~ x / 2 J_i ;0 (i) centered by its mean corresponds 
to Y t N (Ho) —Y N (Ho) since both processes increase (resp. decrease) by iV -1 / 2 
whenever a particle jumps from —1 to (resp. to —1). The limiting covariance 
E[Z s Z t ] corresponds to the formal covariance 



E 



{Y t (H ) - Y (H )}{Y S (H ) - Y (H )} 



Denote by v ^?.\ the measure ^(.) conditioned to have a particle at the origin. 

Remark 2.4. The law of large numbers and the central limit theorem for the em- 
pirical measure and for the current starting from v^*^ follow from the law of large 
numbers and the central limit theorem for the empirical measure and the current 
starting from the measure ^(.) since we may couple both processes in such a way 
that they differ at most at one site at any given time. 

Fix a profile po with the first four derivatives limited, and consider the product 
measure v ^*.y Denote by X t the position at time t > of the particle initially at 
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the origin. A law of large numbers for Xt follows from the hydrodynamic behavior 
of the process: 

Theorem 2.5. Fix t > 0. X t /N converges in P w,» -probability to Ut, the solution 

po(-) 

of 

(d u p)(t,u t ) 

Ut = -, n — • 

P(t,u t ) 

Note that the solution of the previous equation is given by 



dup(t,u) = — ds (d u p)(s,0) . 
Jo 

Theorem 2.6. Assume that po has a bounded fourth derivative. Let Wt = N^ 1 / 2 
(X t — Nut). Under P„n,* , For every k > 1 and every < t\ < ■ ■ ■ < ifc, 

(W^ , . . . , W£?) converges in law to a Gaussian vector (Wt 1 , ■ ■ ■ , Wt k ) with covari- 
ance given by 

p(s,u s )p(t,u t )E[W s W t } = [ dxP Us [B s <x]P Ut [B t <x}xo(x) 

J — OO 
/"CO 

+ / dx P Us [B s > x] P Ut [B t > x] X Q (x) 
Jo 

poo 

dr dxpt~ r {u t ,x)p s - r (u s ,x)xr(x) ■ 



In this formula, P u stands for the probability corresponding to a standard Brownian 
motion starting from u. 

The assumption made on the smoothness of po appears because in the proof of 
Theorem 12. 61 we need a sharp estimate on the difference of the discrete approxima- 
tion of the heat equation (|2.2() and the heat equation l|2.1|l . In section H3 we show 
that there exists a finite constant Co for which [p^ (x) — p(t,x/N)\ < C$tN~ 2 for 
all N > 1, x in Z and t > under the assumption that po has a bounded fourth 
derivative. 

3. NONEQUILIBRIUM FLUCTUATIONS OF THE CURRENT 

Suppose for a moment that the profile po has a compact support. Then, 770 is 
almost surely a configuration with a finite number of particles, and it is easy to see 
that we have a simple formula for the current J_i.o(i): 

J-i,o(t) = J2 r >t( x )- T )o(x) ■ (3-1) 

x>0 

In particular, we can write J_i j o(i) in terms of the fluctuation field: 

-±= |j_ 1>0 (f)-i> n [J-i,o(«)]| = Y t N (H Q )-Y N (H ), 

where H a is the indicator function of the interval [a, 00): 

H a {u) = l{[o,oo)}(u). 

Since the profile has compact support, it is possible to define Y t (Ho) as the limit 
Yt(G n ) for some sequence G n of compact supported function converging to H on 
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compact subsets of R and to prove that Y t N (Hq), defined in a similar way, converges 
to Y t (Ho). 

In the general case, however, when po is an arbitrary profile, neither formula 
makes sense, nor the fluctuation field Y t N (Hq) is well defined. Nevertheless, there is 
a way to calculate the fluctuations of the current by appropriated approximations 
of the function G, as made by Rost and Vares in the equilibrium case. 

Define the sequence {G n : n > 1} of approximating functions of Hq by 

G n {u) = {1 - (u/n)}+l{u> 0} . 



From here we use the next convention: if X is a random variable, we denote by X 
the centered variable X — E v n [X] . 

po(-) 

Proposition 3.1. For every t > 0, 



lim E^N 

n — >oo Po(') 



N- 1 ' 2 J_ 1>0 (t) - Y t N (G n ) + Y Q N (G n ) 



2 



= 



r2 



uniformly in N. 
Proof. Clearly, 

M x . x+1 (t) := J x , x+1 (t) - N 2 I ds{r) s (x)-r) s (x + l)} 

Jo 

is a martingale with quadratic variation given by 

< M XtX+1 > t = N 2 f ds { Vs {x) - r, s {x + l)} 2 



The goal is to express the difference Y t N (G n ) — Y N (G n ) in terms of the martingales 
M XtX +\(i) and to notice that these martingales are orthogonal, since they have no 
common jumps. 
Since 

Jx-i,x(t) - J x , x +i{t) = r)t(x) - Vo( x ) 

for all x in Z d , t > 0, 

Y t N (G n ) - Yq N (G n ) = TV" 1 / 2 G n (x/N){J x - hx (t) - J x , x+ i(t)} . 

A summation by parts and the explicit form of G n permits to rewrite this expression 
as 

nN 

* — ' nN 

x =i 

Representing the currents J x , x +i(t) in terms of the martingales Al x ^ x+ i(t) 1 we ob- 
tain that 

JV-Va j o(t ) _ [Y t N (G n ) - Y N (G n )] 

1 1 1 /"* AT 



We claim that the martingale and the integral term converge to in L 2 (P„w ) . 

po(-) 

In fact, since the martingales are orthogonal, estimating their quadratic variations 
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by tN 2 , an elementary computation shows that 

1 nN 1 

V X — l 



2 f 
< - 

n 



The integral term is more demanding, because in non-equilibrium the two-point 
correlations are not easy to estimate. Expanding the square we have that 



E, 



PO(') 



1 



N Jo 



1 N i 2 

ds — [T} s (0)-rj s (nN)]ds 



2N 



ds / drlK, 



PO(') 



{r)M -rj s (nN))(rj r (0) -rj r (nN)) 



By Lemma f3. 21 the previous expression is less than or equal to Cot 5 ^ 2 n 2 for some 
finite constant Co depending only on p^. This concludes the proof of the proposition. 

□ 

A central limit theorem for the current J_i j o(i) is a consequence of this propo- 
sition. 

Proof of Theorem 12.31 Fix t > and n > 1. By approximating G n in 
L 2 (M)nL 1 (K) by a sequence {H n ^ : k > 1} of smooth functions with compact sup- 
port, recalling Theorem 12. 21 we show that Y t N (G n ) converges in law to a Gaussian 
variable denoted by Yt(G n ). 

By PropositionO {Y t N (G n )-Y N (G n ) : n > 1} is a Cauchy sequence uniformly 
in N. In particular, Y t (G n ) — Y a (G n ) is a Cauchy sequence and converges to a 
Gaussian limit denoted by Y t (H ) — Y n (H n ) . Therefore, by Proposition ^. II J_i i o(t) 
converges in law to Y t (H ) — Y (H ). 

The same argument show that any vector (J-i,o(ii) ; ■ ■ • , J— i,o(ifc)) converges in 
law to (y tl (Ho) — Yo(Ho), . . . ,Y tk (H ) — Y (H )). The covariances can be computed 
since by (|2.3ll 

B[{y t (ffo)-yo(-Ho)}{y.(-Ho)-yo(-Efo)} 

= lim E\{Y t {G n ) - Y (G n )}{Y s (G n ) - Y (G n )} 

n—too L 

= lim { / {(T t G n )(T s G n ) + G 2 n - {T t G n )G n - [T s G n )G n } X o 
n ^°° L Jr 

+ 2 f dr [ (VT t _ r G„) (VT s _ r G n )xr) • 

A long but elementary computation permits to recover the expression presented in 
the statement of the theorem. This concludes the proof. □ 

We conclude this section with some elementary estimates on two points correla- 
tion functions. For < s < t and x ^ y in Z, let 

<f>(t;x,y) = E v n [r) t (x);r) t (y)] , cj>(s,t;x,y) = E v n W] s (x);r] t (y)] . 

po(-) po(-) 

In this formula and below, [/; g] stands for the covariance of / and g with respect 
to [l. 
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Lemma 3.2. There exists a finite constant Co = Co(po) depending only on the 
initial profile po such that 

sup \<j>(t;x,y)\ < ^7^, sup \<f>(s,t;x,y)\ < %-{Vs + J- — } ■ 

x, y ei ^ x, y ez JV 1 \t — s> 

The first statement is a particular case of an estimate proved in jzj- In sake of 
completeness, we present an elementary proof of this lemma. 

Proof. Let Li be the generator of 2 nearest-neighbor symmetric simple exclusion 
processes on Z. An elementary computation shows that cj>(t; x, y) satisfies the dif- 
ference equation 

f (d t <j>){t- x, y) = N 2 (L 2 <j>){t; x, y) -\{\x-y\ = l}N 2 [p N (t, x) - p N (t, y)f , 
\ </>(0;x,y) = 0. 

This equation has an explicit solution which is (negative and) absolutely bounded 

by 

Co(po) / ds¥ x , y [\X s -Y s \ = 1] 
Jo 

for Cq = ||<9po||L- m this formula, (X S ,Y S ) represent the position of the symmetric 
exclusion process speeded up by ./V 2 and starting from {x, y}. A coupling argument 
shows that P x>y [\X S — Y s \ = 1] < ¥^ y [\ X s — Y s \ = 1] where in the second probability 

particles are evolving independently. Since P° y [\X s — Y B \ = 1] < C(siV 2 ) -1 / 2 , the 
first part of the lemma is proved. 

To prove the second statement, recall that we denote by A^v the discrete Lapla- 
cian in Z. <f)(t; y) = 4>{s, t; x, y) satisfies the difference equation 

(d t cP)(t;y) = (A N <p)(t;y) 
(j>(s;y) = (f>(s;x,y) \iy^x, 
4>(s; y) = p N (s, x)[l — p N (s, x)] for y = x. 

This equation has an explicit solution 

<M s ;y) = ^2pt- s (y,z)<i>(s;x,z) + p t - s (y,x)p N (s,x)[i - p N (s,x)} , 

where p s (x, y) stands for the transition probability of a nearest neighbor symmetric 
random walk speeded up by TV 2 . The first part of the lemma together with well 
known estimates on p s permit to conclude. □ 

4. Law of Large Numbers for the Tagged Particle 

In this section we assume the initial measure to be v^f.y the product measure 
^pot ) con ditioned to have a particle at the origin. Keep in mind Remark 12.41 

Fix a positive integer n. The tagged particle is at the right of n at time t if and 
only if the total number of particles in the interval {0, . . . ,n — 1} is less than or 
equal to the current J^i i o(t): 

n-l 

{X t > n} = { J_ 1>0 {t) > J2 Vt(x)} • (4.1) 

x=0 

This equation indicates that a law of large numbers and a central limit theorem for 
the position of the tagged particle are intimately connected to the joint asymptotic 
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behavior of the current and the empirical measure. We prove in this section the 
law of large numbers. 

Denote by \a\ the smallest integer larger than or equal to a. Fix u > and set 
n = [uN] in (|4.1() to obtain that 

{X t >uN} = {7V-ij_ 10 (t) > <7r f,l{[0,u]}> +0(N- 1 )}. (4.2) 

By Theorem 12.11 < , 1{[0, u}} > converges in probability to /„" p(t,w)dw, 
where p is the solution of the heat equation (|2.1ll . 

On the other hand, the law of large numbers for J_i g(t) under ¥ u n is an ele- 

po(-) 

mentary consequence of the central limit theorem proved in the last section and the 
convergence of the expectation of TV -1 J_i.o(i). By the martingale decomposition 
of the current and by Theorem lti.il 

E vN [JV-V-loW] = f d S N[p?(-l)-p?(0)} 

po(-) Jo 

i-t 

d uP (s,0)ds + OiN- 1 ) . 

Hence, TV" 1 J_ 10 (t) converges in probability to — J Q d u p(s, 0)ds. 

In view of 14.2|l and the law of large numbers for the current and the empirical 
measure, 



lim P N „ [N-'X t >u\ = \° [i -{o t duP ^ ° )ds < ■£ P{t ' W)dw 



By symmetry around the origin, a similar statement holds for u < 0. Thus, 
Xf* /N converges to Ut in probability, where Ut is the solution of the implicit equa- 
tion 

p(t,w)dw = — d u p s (0)ds . 
Jo 

5. Central Limit Theorem for the Tagged Particle 

In this section we prove Theorem l2.6l developing the ideas of the previous section. 
Assume first that Ut > and fix a in K. By equation (|4. 1|> . the set {X t > Nut + 
a^/N } is equal to the set in which 

Nu t aVW—l Nu t 

J-i, (t)>J2vt(x) + E Vt(x + Nu t ) - [E^JJ-wW-^pfix)]} , 

x—0 x—1 x=0 

(5.1) 

where pf (x) is the solution of the discrete heat equation (|2.2() . 

We claim that second term on the right hand side of (|5.1|l divided by y/N 
converges to its mean in L 2 . Indeed, by Lemma 13.21 its variance is bounded by 
Coa 2 7V _1 for some finite constant Cq. Since by Theorem lti.il 

aVlV-l 

-= £ p?(x + Nu t ) 

* x— 1 

converges to ap(t,Ut), the second term on the right hand side of (15.1(1 converges in 
probability to ap(t,Ut). 
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An elementary computation based on the definition of ut and on Theorem 16.11 
shows that the third term on the right hand side of l|5.1|l divided by y/W vanishes 
as N | oo. 

Finally, by Proposition ^. II for fixed t, A^ 1 / 2 ! J_i, (i) — Y^,x=o^t( x )} behaves as 
Y t N (G n ) - Y N (G n ) - Y t N (1{[0, u t ]}), as A f oo, n f oo. Repeating the arguments 
presented at the beginning of the proof of Theorem 12.31 we show that this latter 
variable converges in law to a centered Gaussian variable, denoted by Wt, and which 
is formally equal to Y t {H Ut ) — Yq{Hq). 
Up to this point we proved that 

~X t ~ u t N 



lim P^iv 

iV— >oo po(-) 



> a 



= P[W t > ap(t, u t ) 



>N 

provided u t > 0. The same arguments permit to prove the same statement in the 
case Ut = 0, a > 0. By symmetry around the origin, we can recover the other cases: 
u t < and a in M, u t = and a < 0. 

Putting all these facts together, we conclude that for each fixed t, (X t —Nut)/VN 
converges in distribution to the Gaussian Wt/p(t,Ut) = [Yt(H Ut ) — Yo(Ho)]/p(t,Ut). 
The same arguments show that any vector (N~ 1/>2 [X tl — Nu tl ], . . . , N~ 1 / 2 [X tk — 
Nut k ]) converges to the corresponding centered Gaussian vector. 

It remains to compute the covariances, which is long but elementary. 

6. Appendix 

In sake of completeness, we present in this section a result on the approximation 
of the heat equation by solutions of discrete heat equations. 

Fix a profile po : K — > R with a bounded fourth derivative. Let p : M + xl^l 
be the solution of the heat equation with initial profile po : 

d t p(t, x) = d 2 p(t, x) 
p(0,x) = p (x) . 

Recall that we denote by Ajv the discrete Laplacian. For each N G N, define 
Pj (x) as the solution of the system of ordinary differential equations 

" (d/dt)p?(x) = (A N pf)(x) 
p$(x)= Po (x/N) . 

The main result of this section asserts that p N approximates p up to order A~ 2 : 

Theorem 6.1. Assume that po : K — > [0,1] is a function with a bounded fourth 
derivative. There exists a finite constant Cq such that 

C t 



~ N 2 



for all N > 1, t > 0, x € Z. 



An easy way to prove this statement is to introduce a time discrete approximation 
of the heat equation. For each N in N and each S > 0, we define pf' N (k), k in Z, 
I > by the recurrence formula 

/ = pf N (k) + 8N 2 [p 5 { N {k + 1) + pf N (k 1) - 2pf N {k)] 

\ p s > N (k) = P0 (k/N). (6 ' 2) 

We now recall two well known propositions whose combination leads to the proof 
of Theorem 16. II The first one states that the solution of l|6.2|l converges as 8 j to 
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the solution of Ij6.1|l uniformly on compact sets. The second one furnishes a bound 
on the distance between the solution of the discrete equation H6.2|) and the solution 
of the heat equation. 

For a in M, denote by \_a\ the largest integer smaller or equal to a. 

Proposition 6.2. For each N > 1, 

uniformly on compacts o/M + x Z. 

Proposition 6.3. Suppose that 5N 2 < 1/2. Then, there exist a finite constant 
Co = Co(po) such that 

pf N {k)-p(8l,k/N) 

for alll>0,kE Z. 

Clearly, Theorem 16.11 is an immediate consequence of Propositions 16.21 and 16.31 
Proposition l6.2l is a consequence of Proposition ^ . 3l and the Cauchy-Peano existance 
theorem for ordinary differential equations. Proposition 16 . 31 is a standard result on 
numerical analysis (see 9 for instance). 
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